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Abstract 

The propagator of linear molecules whose constituents interact 
through oscillator potentials can be obtained in a closed form for N 
atoms as long as < 4. We compute the propagator for arbitrary 
N in the approximation ^ 1. Taking advantage of this result it 
is possible to analyze the limit in which the molecule has an infinite 
number of constituents with infinitesimal length of sepparation, cor- 
responding to the quantization of a string, elastic rod or the second 
quantization of a Klein Gordon particle. The evolution of some spe- 
cific initial conditions is also studied, namely the time development of 
states with minimal dispersion and the effect of sudden perturbations 
on the vacuum of the scalar field theory. 
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1 Introduction 

Transient phenomena in quantum mechanics have been of interest since the 
developments in [1] a long time ago and others which are more recent [2], 
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[3]. The treatment of this subject demands a dynamical study of different 
quantum systems corresponding to specific problems. In view of the simple 
theoretical results obtained for models involving oscillators in [4], there arises 
naturally the question of their generalization to systems with infinite degrees 
of freedom. 

In many textbooks [5] we find treatments of continuous media as the 
limit of interacting discrete systems in the context of classical mechanics. 
However, it is possible to exploit this analogy in the quantum picture through 
the Feynman formulation of the propagator [6], connecting the latter with 
classical lagrangians. In this fashion one may study discretized versions of 
second order field equations and quantize them before taking the continuum 
limit. The inverse order of this limits makes it possible to perform certain 
calculations, while when taken in the usual order gives rise to path integrals 
of field configurations. The latter may be difficult to handle or may even lack 
of an appropriate mathematical definition. 

Our prototype discrete system will be the linear molecule and we proceed 
to compute its propagator in section 2. In section 3 the continuum limit of 
our system is taken calculating thereby the propagator for strings and second 
quantized scalar fields. In sections 4 and 5 we apply our results to certain 
initial conditions, including sudden external interactions. 



2 Propagator of the linear molecule 

We start with a lagrangian in 1 + 1 dimensions of a chain of N particles of 
equal masses m connected by equal springs of constant k — muj'^ 



>C=?E(i?-'^'(^.-^mr)-E>C. (1) 
^ 1=1 i=i 

It is costumary to find the orthogonal transformation which sepparates the 
problem in N independent lagrangians (normal modes) and we will do so in 
the following. Letting (x)^ = Xi and 
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our lagrangian becomes 
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L =^(x2-a;2x^Vx) (3) 

and it is left to diagonalize V. It can be shown (see Appendix) that the 
eigenvalues of V are given by 



A„ = 4sin^(^ ^^^_^J +Q(l/Ar), < n < TV - 1 (4) 

which is particularly useful for large N . We can see that ~ 10 is already 
a good approximation as can be checked numerically. We resort to this type 
of formula since it allows to manipulate our expressions simbolically though 
numerical methods to compute exist. Notice that in the limit iV — > oo 
the A's (frequencies of the normal modes) populate densely the interval [0, 4] 
distributed by a trigonometric law. It is straightforward then to write 



V^Dy) 



(5) 



where y = Ox, O being the orthogonal transformation corresponding to the 
diagonalization of V into D = diag (Aq, An-i). The propagator is written 
as [7] 



mouX 



1/2 



n=0 \ 2mhsm[LVAn t 



X exp 



2^ sin 



cos {uj\]lh){yl + y'l)-2yny'^) 



(6) 
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and this expression is exact as long as the A's contain corrections 0{1/N), 
but these do not alter the functional form of (6). Notice that the factor 
associated to n = corresponds to Aq = and can be written as 



-f^o(yo,y'o;^) = 



exp 



im . 



y'oY 



(7) 



27riht 

which is the propagator of the free particle, finding thereby that yo is the 
coordinate of the center of mass. The meaningful degrees of freedom are 
therefore indicated by y, but we can at any moment employ the inverse 
transformation x = O^y since O can be written explicitly (see Appendix) . 



3 The continuum limit of the hnear molecule 
propagator 

Here we will consider the continuum limit as that indicated in [5] or [8] . First 
of all let us recall that the lagrangian (1) can be multiplied by a parameter 
a with the dimensions of length so that 



^ /m.2 rmj'^a Xi- Xi+i^\ 

^ r dAd,<pr - ^{d,<^r (9) 

iV-*oo, a-+0 Jo 2 Z 

as long as Xi(t) a = — > c?^, ua — > fl. The constant L is 

interpreted as the length of the resulting string and fl is such that /xO^ = Y, 
ie. Young's module. 

For the propagator limit we use Feynman's prescription 



r„ C" (t /„' r '<'''*^' - ^('^*^') 

which is the propagator of a string, elastic rod or a massless scalar relativistic 
particle when the classical length L — > oo and the identification Q <-> c takes 
place. 
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Though we are quite famihar with this form of the limit . in normal coor- 
dinates y some subtleties arise. First of all, we notice that the ^ oo limit 
of lagrangian (5) is equal to the limit of the following lagrangian 

E«(f^|-^(A./a^ (11) 
where we can observe that the eigenvalues are now \j/a^. Since 

we see that these eigenvalues are distributed discretely for finite j by a square 
law, just as the energy spectrum of a string with given initial and final 
conditions at a fixed time. Therefore, the limit of (11) is not an integral, 
but a series and it will make sense if we set ^/ayj — > 77^, ie. if fields are 
rescaled. To this effect we should also rescale the orthogonal operator O in 
the following way 

^O,, ^ 0,(0 (13) 
so that the transformations between fields become 



JO ■ 

fo,{i)Oi{i) = E 0,(00,(0 = - n (14) 

JO ■ 

where a'^Sai — > (5(^ — {'). With all these considerations, lagrangian (5) 
posseses the limit 



>C - E - ^{j^/L)W, (15) 
i=o ^ ^ 

whose discrete form stems from the discreteness of the eigenvalues A, as we 
have proven. 

The propagator related to lagrangian (15) can be written directly or it 
can be computed as the limit N ^ 00 oi propagator (6) but taking into 
account that the prefactor is proportional to m^/^ ~ oJ^I'^yJ^I'^ . Thus, for the 
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fields r]j the propagator is multiplied by the jacobian a cancelling the 
scale factor emerging from the mass. The resulting expression is 



with 



2% ^ Lsin(Q_77rt/L) 



m ^ n 



+ cos(Oj7rVL) - 277,-77^.] ) (16) 



Q 2mhLsm{ijjjT:t/L) 



(17) 



The definition for / deserves some comments. The convergence of this 
infinite product can be spoiled by the presence of and leading to f — 
cx), 1 or for different values of these parameters. Nevertheless the propagator 
itself obeys a normalization condition even when the prcfactor is a divergent 
or a vanisliing quantity. This stems from the fact that K is the representation 
of a unitary operator and preserves norm. The limit of propagators should 
be taken under the integral / D(f) [•] or its discrete version from which it is 

easy to see that a rescaling of fields 4> = \J V = \f^V g^^s rid of the 
problem. Using dimensionless fields t) is the price to be paid. With all this 
the propagator becomes 



oo 



-f^string {v,V';'t) = /(^)exp 



^^§Lsin(0j7rt/L) 



{r]] + Vj') cos{nj7rt/L) 



with 



m ^ n 



jn/L 



j=0 



i sm{ujj7rt/L) 



det 



i{V + ie) (sin(V^ + ie))" 



(18) 



(19) 



V = Voo and e arbitrarily small. Let us now write the propagator of a second- 
quantised Klein Gordon particle of mass M. First of all we need to add a 

1). Since this term is invariant under 



mass term ""^^'^^^ x • x to lagrangian 



2h 



the orthogonal transformation with associated matrix O, it is sufficient to 
replace u in (6) by uij = \ oj"^ + ^® have indicated, we must set 
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Q — c. Once this is done we need to extend the molecule to the whole real line 
by taking L oo. We identify the continuous wave number as jn / L k, 
giving a differential n/L ^ dk. The fields fj must be rcscaled again in 
order to have the limit \/Lf]j — > vif^)- The orthogonal transformation obeys 
VLOj{$,) 0{^,k) in order to have field transformations in full integral 
form, i.e. 

dcoik,om, dkoickuk) 

-oo J — oo 

/oo roc 

d^O{k, oo{k', - S{k -k'), / dkOi^, k)0{^', k) - 5(e - 0(20) 

-oo J —oo 

where L5jji — > 5{k — k'). Meanwhile, the lagrangian (15) with a mass term 
has the following limit when L — > oo 

/oo 
dkm)\^-{E,/hf\r^{k)\^ (21) 
-oo 

where we have computed 

h\]'^Uj ' \ffF^]^^^M^ = Ek (22) 

From the equations of motion induced by this lagrangian we see that 
spatial derivatives become factors of k, inferring thus that 0{k,^) is the 
Fourier kernel. Finally the sought propagator is given by 



g{t) exp 
with 



/: 



iEk 



cos 



h 



Kkg [r], r]', t] 
{\r^{k)\'+W{k)\')-2r^{k)r^\k) 



(23) 



git) ^ n 



E, 



(24) 



The prefactor is given again in terms of an infinite product and it is left 
left indicated without discussing its convergence. At any rate, it seems more 
suitable to work with the molecule model and its propagator (6) or its mas- 
sive version in order to study dynamical effects on strings or particle fields. 
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Any meaningful time dependent quantity of these systems will be treated 
discretely and limits will be computed at the end of calculations. 



4 Minimally localized molecule and its con- 
tinuum limits 

At this point we are ready to apply propagator (6) to some initial conditions. 
Consider the linear molecule of atoms to be in a state which is minimally 
dispersed, i.e. a gaussian distribution in space for each atom: 



/ 1 \^/^ 



N-1 



4a2 



1=0 



(25) 



where the width a is the same for all components. This wave function has 
the property 



^(x',0) = 7A(Ox',0)=^(y',0). (26) 

and it describes a localization in phase space with minimal uncertainty. In 
the continuum limit A^ oo, initial condition (25) becomes 



where cr^a s^. Similar remarks to those following (17) hold here. Despite 
the limit of the factor preceding the exponential may be vanishing or diver- 
gent, functional (27) is normahzed and for all t we have / D(t)\ip [0, t] p = 1. 
Thus, functional (27) is a probability amplitude depending on field configu- 
rations. In order to obtain a normalized functional in normal coordinates we 
need though to use the jacobian a~^/^ to write the initial condition as 



so that 



Now we apply (6) to (26) by integrating over y' to obtain 



1 



'N-l I 

5 V2^^|W, 



1/4N 



exp 



N-l 



E 



4a|(t) 



with 



aj{t) = a 



\ 



1 + 



and A(t) a phase which is irrelevant to the probability density: 



wy,^)r-ln(^) I exp 



(29) 



(30) 



- - 1^ sin^ (uXyh) (31) 



(32) 



We observe therefore an oscillatory behaviour of the density widths and am- 

1 / 2 

plitudes depending on the corresponding frequencies uXj . The center of 
mass (Ao = 0) distribution exhibits disispation as expected. In the contin- 
uum limit (31) yields 



N—*oo 



1 + 



- 1 sin^ {^tjn/L) (33) 



whose maximum value decreases with j. For the massive Klein Gordon par- 
ticle we find a different hmit. The discrete case yields 



aj{t) = a 



1 + 



-l]sm^{u;jXyh) 



(34) 



Taking ^ oo, L — > oo, initial condition (28) becomes 
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lim , 



Af/4- 



1 



(35) 



after rescaling widths s = VW*" '^^^ wave functional for fields r]{k) at time 
t is now 



V5 [77, t] 



1 



lim , 



—00 



(36) 



with 



= 



1 



'2^ 



sm 



n 



(37) 



This quantity is again oscillatory but with a frequency given by the relativis- 
tic formula of kinetic energy. The most probable field configuration remains 
7^ = 0, but the field localization width can spread to a factor {-s/^Tic) / {EkS^). 

To illustrate the significance of this result let us compute the average 
number of quanta for each mode using states -0(7, t)- In the discrete case we 
have 



d 



oscillator (^-i^, wA/ J - 1/2 ) (38) 



(39) 



Taking the continuum limit N ^ oo the average number results in 



with Sj{t) given by (33). For the massive relativistic particle we have to 
replace the frequencies by ujj and take L — > oo. The expectation values 
computed by integrating with respect to fields r){k) are now multiplied by 
the squared jacobian 1/L ^ dk/ir, thus we get 



f4W + ^(4W)-'-i 



rik 

- ^ iy,it)dk (41) 

TT 
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where we have obtained a number density Vkit) in momentum space as a 
function of time. 




Figure 1: Average number density of particles as a function of t for k — 
(Solid hne) and A; = ^ (Dashed hne). 

An oscillatory density number is not a surprising result when viewed as 
the evolution of harmonic oscillator number of quanta. However, it is inter- 
esting to notice that such density number corresponds to a second-quantised 
free field which has beeen set initially in the most localized configuration 
in phase space (gaussian distributions). 

5 The string under a sudden electrostatic field 
and the relativistic particle with a source 

Now we may think of another example which involves a scalar field and 
a sudden interaction with a source which is equivalent to a charged string 
between the plates of a condensator in the discrete case. Consider a system 
whose lagrangian is given by 
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L' = L +e{t)£ -x 



(42) 



where C is given by (1), 6{t) is the step function and S denotes some elec- 
trostatic field multiplied by the charge of the ith particle in the molecule, 
i.e.Si = QiEi (remember that this is a one dimensional example and vectors 
here stand for the N components of our system). If we set ^ p{^) as the 
hnear charge density and ^ p{^)E{^) = £{^), then the limit of (42) is 
well known 

^' ~^ r^^ef^w(e,t))'-^(5«0(e,t)r + ^(tK(O0(e,t)') (43) 

N-^oo Jo \Z Z J 

and leads to equations of motion [8] 



d^cf>-n^dlcl>=JiO, t>0 (44) 

where we have defined a current J(^) = (0- But before taking the con- 
tinuum limit, let us assume that for i < the discrete system is in its ground 
state, i.e. 



N-l 



1/4 



^(y,o)= n 



i=0 



exp 



1/2 



|0,...,0) (45) 



given, of course, in normal coordinates. The source term in (42) can be 
written as 



£-x = £'-y, £' = OS (46) 

With this orthogonal transformation of the external field and for positive 
times we can immediately write 
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/:'(y,y) = ^(z,z) + E 



i=0 



(47) 
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C I 

where Zi = yi — ^ji^. ■ In the N ^ oo hmit we must set £'j/ ^Ja — > Fj so that 



Jo 



(48) 



The propagator related to lagrangian (47) has the form (6) but evaluated in 
z and with an additional (but irrelevant) phase factor containing the field 

energy. 

It is quite easy to apply such propagator to the vacuum state (45) and 
we write the resulting probability density computed in [4] for each mode 



1/2 



1/2 



exp 



.1/2 



(49) 



where we can see a time dependent average of the distribution given by 



aj{t) 



The continuum limit yields 



(50) 



aj{t)/y/a 



FjSm^{ntjTx/{2L)) _ 



(51) 



For the Klein Gordon Field we again rescale fields and set, therefore, VLFj — > 
F{k) so that F{k) can be the Fourier transform of S {^) when L — >■ oo. The 
average field configuration (51) in the massless case yields 



llJ,nL 



2h "^-^'^ 



— > , ^ , sin^ (ct/c) 

a{Kt) = I{k) J ^ 



(52) 



where I{k) is the Fourier transform of J{C)- As before, the massive case is 
obtained by replacing the frequencies and the result is 



a{k, t) 



h^I{k) sin' (1^) 



(53) 
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The average number taken between states (49) yields, in the discrete case, a 
number density in momentum space after taking the hmits N oo, L oo. 
Proceeding analogously to the derivation of (41) we arrive at 



1 _^ h^\I(k)\^ sm^(^) 

jm ^' ^'''U k^u,it)dk (54) 

In this result wc observe that the Fourier transform of the current {I{k)) 
determines whether the number density u increases with energy or not, de- 
pending on the explicit form of J(0- Another unexpected feature is that even 
when the source is constant in time (say t > 0), the number oscillates with 
a frequency depending on the relativistic energy of the produced particles. 



6 Conclusions 

Discrete systems are useful to study quantum dynamical problems whose con- 
tinuum limits may be difficult to formulate in quantum field theory. We have 
outlined a method by which we can obtain propagators (6,18,23), function- 
al (30,32, 49) and physical quantities (40,41,54) by developing some simple 
examples: molecules, strings (or elastic rods) and Klein Gordon particles in 
second quantization. The method allows interactions in a non perturbative 
regime, but their form is kept simple so that lagrangians remain quadratic 
in their discrete version. 



Appendix 

Now we proceed to derive expression (4). Consider the matrix V^v as in (2) 
and let 



M 



N 



/ 1 -1 
-12-10 
0-12-1 








-1 




\ 





2 -1 
-1 2 J 



(55) 



NxN 
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Define also the polynomials 



0iv(A) = |V^- AIatI, x^(A) = |Mjv- AIjvI (56) 
It can be easily verified that the following recursion relations hold 



/ M^) \_( I -A W 07v-i(A) A . . 

[xn{X) J [l 1-X )[xN-i{X) J ^^'^ 

The matrix appearing in the RHS of (57) can be diagonahzed and then 
the recursion is solved by computing the powers of that matrix. Let a — 
|(A + ^Jx{X--4)) and choose as initial conditions 02, X2 which are given by 



02(A) = A' - 2A, X2(A) = A' - 3A + 1 (58) 

Following the procedure indicated in the last paragraph by finding the ap- 
propriate similarity transformation we can write (57) as 



(1 - a)^-2 A 

(1 - A)A^-2 ) X 

^ ( -:t ) ( tl^', ) 

from which the first component is used to write the secular equation 0Ar(A) = 
0. We notice that the reality of A implies a = e*" and therefore A = 
4cos^ {a/ 2). After some algebraic steps we find that the secular equation 
is equivalent to 

(02(A) - 2x2(A)) sin {{N - 2)a) cos (a/2) - 02(A) cos ((iV - 2)a) sin (a/2) = 0(60) 
Now we can estimate the solutions of (60) by making the Ansatz 



a — 



N -1-n 
N-1 



-TT 
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(61) 



for finite n. Replacing back in (60) and neglecting terms 0{1/N) we see that 
the secular equation is solved. Therefore 

A = 4sin2(^^^^)+0(l/iV) (62) 

We can also write the equations for the eigenvectors in recurrence form and 
these can be solved in terms of the eigenvalues given above. 
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